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Abstract 

The Hamiltonian of a moving atom in electromagnetic fields includes velocity-dependent 
terms. We show that the leading velocity dependence emerges systematically in the non- 
relativistic limit from a scheme firmly based on the relativistic invariance of the energy- 
momentum stress tensor of the coupled matter-fields system. We then extend the Hamiltonian 
to the many-body situation suitable for describing a Bose-Einstein condensate (BEC). From 
first principles, we use the equation of motion for the condensate wavefunction to obtain an 
extended version of the Gross-Pitaevskii (GP) equation and an equation for the internal states 
of the atoms. It is shown that laser fields modify the GP equation by inclusion of convective 
terms involving a Rontgen interaction plus a term coupling the centre of mass momentum to 
the Poynting vector. We also obtain the modified Maxwell equations for the electromagnetic 
fields coupled to the BEC involving the average velocity of the atoms. 
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1. Introduction 

In a recent letter, Leonhardt and Piwnicki [1] examined magnetic effects associated with the 
motion of a Bose-Einstein condensate, regarded as a quantum dielectric. Earlier, it has been 
pointed out that a neutral atom bearing an electric dipole moment moving in an external 
magnetic field can accumulate a quantum phase [2-6] and a moving dipole may, under 
suitable conditions, exhibit a detectable Aharanov-Bohm phase shift [7]. These and other 
effects associated with atomic motion have received considerable attention recently, most 
notably in the areas of atom optics [8,9] and laser cooling and trapping [10]. The theoretical 
framework describing relevant phenomena involving moving atoms relies on a seemingly 
straightforward extension of non-relativistic quantum optics by including the translational 
motion of the atomic centre of mass. Clearly, it is this feature of the theory which ensures that 
the exchange of energy and linear momentum is properly taken into account in any 
interaction between the laser light and the atomic translational and internal degrees of 
freedom. 

The need to incorporate the centre of mass motion in quantum optics theory had necessitated 
a re-appraisal of the corresponding quantum electro-dynamical theory where investigations 
sought to elucidate how the division of the motion into centre of mass and internal motions is 
affected by the presence of the interaction with electromagnetic fields [11-14]. One of the 
main outcomes of these investigations was the emphasis on the role of the Rontgen 
interaction energy term which couples the electric dipole moment to an effective electric field 
involving the centre of mass velocity and the magnetic field. Although the existence of the 
Rontgen interaction in quantum electrodynamics was pointed out earlier, its practical 
implications and its significance for the consistency of any treatment incorporating 
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translational effects in the coupling of electric dipoles to electromagnetic fields have only 
relatively recently been explored [1-7,15-17]. 

Our first goal here is to incorporate the motion of the centre of mass in a systematic manner, 
aiming to derive a multipolar Hamiltonian which exhibits its velocity dependence. This is 
safely done using a relativistic starting point. Having obtained the velocity dependent 
multipolar Hamiltonian, we take the non-relativistic limit and incorporate the many body 
effects appropriate for describing a Bose-Einstein condensate (BEC). We must also add the 
S-wave scattering potential to the Hamiltonian to take into account the weak hard sphere 
interactions which occur in a dilute BEC [18,19]. In a similar manner to other works [18], we 
find the equation of motion for the condensate wavefunction assuming that the majority of 
the bosons are in their ground state. This leads to a modified Gross-Pitaevskii equation which 
now incorporates velocity dependent terms in the coupling of the atoms to the 
electromagnetic fields. As stated elsewhere [20], we must be careful how the polarization 
terms are included and divide the Coulomb interaction terms into intra- and inter-atom 
interactions. After taking the dipole approximation [21] and assuming that the electron 
position vector can be averaged in the interaction term, we may uncouple the no n- linear 
Schrodinger equation into an equation for the gross motion of the atoms and another for the 
internal degrees of freedom. 

In [1] it was shown how the electromagnetic fields can be modified in the presence of the 
BEC by including the average velocity of the atoms. Here we show how Maxwell's equations 
can be obtained as Heisenberg equations using the Hamiltonian we have derived and show 
that the sign of the velocity in [1] is different an observation which can be justified further by 
inspection of the Rontgen interaction. We conclude by breifly discussing the results obtained. 
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2. Transformation of the Hamiltonian 

We start by writing the fully relativistic Lagrangian of the system in a frame of reference 
moving with the atom including the electromagnetic field and its interaction with the Dirac 
field for a bound electron. A contribution due to the massive nucleus is included, with the 
nucleus taken to coincide with the centre of mass, providing the source of the binding 
Coulomb field and is assumed to be moving relative to the laboratory observer. From the 
Lagrangian in the centre of mass frame it is straightforward to derive the stress energy- 
momentum tensor T pv , with the T 00 tensor component identified as the Hamiltonian of the 
system in the centre of mass frame. Under a Lorentz transformation from the centre of mass 
frame to the laboratory frame the stress energy-momentum tensor transforms like a second 
rank four-tensor 

T^ v =Af x a A v p T a P, (1) 

where A^, is the Lorentz transformation matrix with an arbitrary relative velocity vector V. 

The required Hamiltonian in the laboratory frame emerges as I 00 and can be written as a sum 
of a contribution due to the electromagnetic field, a contribution due to the Dirac field 
together with the centre of mass term and, finally, a contribution describing interactions. The 
velocity dependence is manifest in various terms of this Hamiltonian. 

The system Hamiltonian in the laboratory frame is next subjected to two sets of 
transformations. First, an exact mutipolar gauge transformation is applied and, second, the 
non-relativistic limit is taken. In the process of dividing the resulting Hamiltonian into 
unperturbed and interaction terms we retain terms to first order in the velocity. This velocity 
is identified with the momentum operator of the centre of mass. 
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The result can be written as 



H=H f +H a +H int . (2) 

The first term is the usual electromagnetic Hamiltonian plus a term which represents the 
kinetic energy arising from the transport of the field Poynting vector at velocity V, 

H f =^|[£ l2 (r) + c 2 5(r)pr-^|p.k(r)xB(r)t 3 r. (3) 

Here E 1 is the transverse part of the electric field and P is the centre of mass momentum of 
the atom which is given by V / M. The second term contains the gross kinetic energy of the 
atom, the kinetic energy of the electron and the polarization term containing the interaction 
between them in the Power-Zienau-Woolley representation [21], 



2 2 

H * = ^7 + J h + ^~ \p 2 (r,*,\)d 3 r, (4) 
2M 2m 2e J 



' (r'R'l) = e|l(5(r - R - ti)dA, (5) 


p is the momentum of the electron which is at position R + 1 with the centre of mass of the 
atom (assumed to be at the position of the proton) at R. The last term in Eq. (2) contains the 
interaction between the atom and the electromagnetic fields, 

Hint =~\' 1 (r,R,l)E ± (r>/ 3 r + ^J'p{" 1 (r,R,l)xB(r)^ 3 r. (6) 

In the above we have disregarded all magnetic multipolar interactions as well as spin 
interactions but the full electric multipoles as in Eq. (5) are retained at this stage. 
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3. The condensate wavefunction 

We now transform the Hamiltonian in Eq. (2) to the case where there are many atoms, 
suitable for the case of a Bose-Einstein condensate (BEC). The boson field operator 
describing the internal as well as the gross motion of the ensemble can be written as a sum 
over the product states as follows, 



^ t (R,l)= ^y/^Rk^l)^, 

nl,n2 

*(r,i)= £v^(Rk n2 (i>W 



(7) 



nl,n2 



where y/ n i(r) are the state functions associated with the centre of mass of the atom, Xn2(\) are 

associated with the internal states of the atom and the operators n2 and a nl n2 are the 

usual boson operators satisfying the usual commutation relations. Thus, instead of Eq. (2) we 
have the many-body field theoretic Hamiltonian 



H -\\ 



^(RT) 



D 2 2 
P p 

+ — 

2M 2m 



*(R,1)+ 

^ r |[ v i ,t (R,l)P v i , (R,l)][' 1 (r)xB(r)-e E 1 (r)xB(r)J^ 3 r + 

i||^ t (R,l)^ t (R,iy (R-R , >i , (R , ,i , ) x l'(R^ , V 3 ^ 3 ^+ 
^{R,\y tra JR)V{R\Y) 



d 3 Rd 3 l + 



2e 



d 3 r, 



(8) 



where we have introduced the polarization field operator (r) in the form 



{r)=^e^Yi lt {R,l)5{r-R-M)i t {R,l)dXd 3 Rd : 



(9) 
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We have added two terms to this Hamiltonian. The term involving Vo(R-R') i s tne 
scattering term due to hard- sphere interactions [18,19]. This is usually taken to be in the S- 
wave approximation [18,19] which is 



iTth 2 a 



Vq ( R _ R ') = ^± d (R - R') = v 8 (R - R'), 



M 



(10) 



where a is the scattering length. The second term we have added is due to the trapping 
potential V tra p(R) which is usually taken to be parabolic. 

We now derive the equation of motion for ^(Rj) using the Heisenberg picture 
-ifW = [#, V I / ]. Dropping all velocity dependent terms and terms containing the polarization, 
we have , using ^ = -iE/h^P where E is the total energy of the state, 



+ f" + Vtrav (R) + V V f (R, 1)F(R, 1) 



2M 2m tmp 



^(R,l) = ^(R,l), 



(ID 



which, apart from the implicit dependence on the internal states ^(1) and the kinetic energy 
term for the electrons, is reminiscent of the Gross-Pitaevskii (GP) equation prior to assuming 
that all particles are in the ground state. 

We already see that Eq. (11) requires at least the Coulomb interaction between the proton and 
electron of the same atom to account for the correct internal states. This Coulomb term comes 



from the longitudinal part of the squared polarization term, 



if ^ 



r, in the 



Hamiltonian. The transverse part, — — f 1 (r)j 3 r, amounts to a self-energy of the atom 

2e J 

which we can ignore. We must, however, add the Coulomb interaction between the proton 
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and electron to the full Hamiltonian, Eq. (8), and subtract it. After some lengthy algebra, we 
obtain the new form for the Schrondinger equation, which on restricting to the electron dipole 

approximation where ' (r)= e||l v I /t (R,l)5(r-R) v i , (R,l)j 3 Rd^l , becomes 



D 2 2 2 
P p e e 
- + - + 



Jd.!'^ 1 (r, r)p(R, r) - y(r, i'Jv 1 (R, l')J/ 3 /' 



+ 



2M 2m 4ne l e 
y fra/? (R)+y ^ t (R,l)^(R,l)-d.E 1 (R)+[ixB(R)-£ E 1 (R)xB(R)]^ r 



^(R,l)=£ v i / (R,l). 



(12) 



where we have put d = el. Next we make the assumption that the number of atoms in the 
condensate is very large and that they all occupy the ground state without any fluctuations. 

Thus, the ground state operators reduce to c-numbers, ciq q ~«qq ~ y[N , with all other 

operators a Ml? „ 2 n (nj,n 2 > 0) giving zero. Even so, in Eq. (12) the internal and gross 

motions are still coupled. To uncouple the two motions, we assume that 1 can be averaged in 
the interaction term and the dipole moment. We then get the total energy as E = fi + E mi , 

where fi is the chemical potential of the atoms appearing in a modified Gross-Pitaevskii 
equation, 



2M 



■ + 



Jd.l'| Zo (lf d 3 l'+V 
dE 1 (R)+ [l x B(R)- EqE 1 (R)x B(R)] 



e 

£0 



+ V trap {R)- 



trap 

P 

M 



Vo(R)=A*Vo(R), (13) 



and E- mt is the internal energy of the atoms which appears in the equation describing the 
internal degrees of freedom, 
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2m Ane^l 



*o(l)=£int*o(l). 



(14) 



If we assume that there is a maximum correlation between the atoms, Eq. (13) becomes 



2M 



+ 



( ? 



+ V 



+v trap {R)- 



d.E 1 (R)+ ^lxB(R)-e E 1 (R)xB(R)] 



P 

M 



VA (R)=M^o( R ) 



(15) 



Thus, the dipole-dipole interaction term and the last three terms in the bracket are new terms 
in the non-linear Schrodinger equation. The dipole-dipole term, shown here to arise 
systematically from the theory., has previously been added phenomeno logically [20,22]. 
However, to the authors' knowledge the P-dependent interaction terms involving E and B 
fields have not been derived rigorously before. 

4. The modified Maxwell equations 

We now turn our attention to the electromagnetic fields and the Maxwell equations associated 
with them. To this end we use the Heisenberg picture with the operators B and D = £oE ± and 
the commutation relation [23], 



^ Et(r),Bj(r')]= -ihe ijk -Kd{v-v'\ 



(16) 



where % is 1 if ijk is cyclic, if any of the indices are the same or -1 otherwise. We use this 
commutation relation because the fields are those in the laboratory frame in the transformed 
Hamiltonian and they have not coupled to the polarization. Thus we obtain the following 
Maxwell equations 
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VxH = 



3D 



VxE = - 



3B 

3f * 



(17) 



where the fields H and D contain polarization and velocity dependence as follows, 



H(r)= J u B(r)-e ^VvxE ± (r), 



(18) 



D(r) = £ E(r)+~ 1 (r)-e /VvxB(r), 



(19) 



with the average velocity of the atoms given by 



1 



J v ^(R)Py'o(R^ 3 R 



(20) 



v = — 



M 



In the electric dipole approximation, and assuming that all the atoms occupy the ground state, 
the polarization, Eq. (9),is given by 



Eqs. (17-21) are similar to the Maxwell equations derived in [1] except for the sign of the 
average velocity vector. It should be noted that the starting point in [1] os a model 
Lagrangian and the equations derived from it have the Rontgen interaction with a different 
sign from that published elsewhere [11,13-15]. Here, the equations arise from a theory based 
on a well-know Lagrangian and the associated Hamiltonian which we then systematically 
reduced to yield a theory for a BEC in the presence of external fields. 

5. Conclusions 

We have rigorously derived the non-linear Schrodinger equation governing the ground state 
wavefunction of a BEC including the dipole-dipole interactions and the internal degrees of 
freedom of the atoms. The interaction of the atoms with external electromagnetic fields is 



(r)=/Vd|v (rf. 



(21) 
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derived here including a BEC for the first time, appearing in direct and velocity dependent 
terms. We have also obtained the associated Maxwell equations for the electromagnetic fields 
and shown that these equations include terms dependent on the average velocity of the atoms 
and their polarization. The influence of the new interaction terms on the properties of a BEC 
in external fields has yet to be investigated, but it appears that they become significant when 
vortices are present [1]. This will not be pursued any further here. 
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